Phase information revealed by interferences in the ionization of rotational wave 

packets 
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Time-resolved photoelectron spectra are proposed for the measurement of classical information 
recorded in the quantum phases of a molecular rotational wave packet. Taking Li2 as a prototypical 
system, we show that an interference arises from the electron-nuclei entanglement induced by the 
molecular anisotropy. This phenomenon is used to transfer the information that was stored initially 
in the nuclear rotational degree of freedom into the electronic degree of freedom. 
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I. INTRODUCTION 



In parallel with quantum measurement and quan- 
tum entanglement, quantum interference is one of the 
usual ingredients of quantum communication and com- 
puting The overall speedup expected from quan- 
tum computers originates in part from the fact that the 
input state of the computation can be chosen as a co- 
herent superposition of all possible classical inputs. A 
sequence of unitary operations is applied to this initial 
state in order to perform the computation. This ap- 
proach processes the data in a massively parallel way, 
and the result of the computation usually depends on 
the interference between the various paths followed by 
all possible initial states. The real power of quantum 
interference is revealed when bringing together all differ- 
ent components of the wave function in a single step '5| . 
Even if not sufficient for reliable and efficient quantum 
computation 0, Q , the buildup of a quantum interfer- 
ometer is of course necessary for the construction of a 
quantum computer. In this respect various physical im- 
plementations have been proposed to test the principles 
of quantum computing, including nuclear magnetic res- 
onance 13, trapped ions cavity quantum electrody- 
namics 6], Josephson junctions 7] and neutral atoms [g. 

Matter- wave interferometry using molecular systems is 
an active field of research with potential applications 
for quantum information processing. In this article, we 
show that in a molecular system quantum interferences 
can be used in combination with electron-nuclear entan- 
glement to extract classical information initially stored in 
a rotational wave packet. This approach is inspired by a 
set of beautiful ex peri ments performed on Rydberg atom 
data registers [lOl [Tj , and on Li2 in the context of fem- 
tochemistry |l2j| and for coherent control [T^ [T3 | . This 
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last system has also been used for the implementation of 
various simple quantum algorithms [l5l |. 



II. THE MOLECULAR SYSTEM AND THE 
LASER FIELDS 

The present theoretical study is based on the usual 
three-step photoionization scheme used in most recent 
experiments 12, 13], as shown in Figure ^ 

A cw linearly polarized laser pulse first excites a sin- 
gle transition from the ground electronic state X {^^g) 
of Li2 to a pure rovibrational {va, Na) level on the first 
excited A (^Su ) electronic potential curve. From this 
unique launch state an ultrafast laser pulse of linear po- 
larization excites a superposition of rovibrational levels 
{ve,Ne = Na±1) on the E elect ronic potential. 

We assume here that this pump pulse has been sent 
through a dispersion-free pulse shaper. This kind of tech- 
nique has already been used with Li2 [13| . A phase mask 
is then used to control accurately the phases of the var- 
ious components of the nuclear wave packet created on 
the E electronic state potential. 

Each possible value of the projection AI e [—Nx : Nx] 
of the initial state rotational quantum number on the 
electric field polarization axis is equiprobable. We there- 
fore perform a separate calculation for each initial value 
of M, and the total photoelectron signal is calculated 
by averaging our results over this initial state quantum 

Note that the three excitation 



number (section IIII D|) 
steps shown in Figure ^ are performed with the same 
electric field linear polarization, and the value of M is 
thus identical in all electronic states of Li2 . 

In analogy with the study performed on Rydberg 
atom data registers |0, |^ , we associate the phases im- 
printed in the coherent superposition of rovibrational lev- 
els {ve, Ne) with information of classical nature. For ex- 
ample, if a single vibrational level w_b is populated, the 
phase shift Aip between the two rotational components 
Ne = Na + 1 and Ne — Na — 1 is fixed at the value 
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FIG. 1: (Color online) Three step "preparation-pump-probe" 
scheme proposed to store classical information in the phases 
of a rotational wave packet in the E (^S^) electronic state 
of the Li2 molecule, and to read it in the photoelectron spec- 
trum recorded in the last step which consists of a picosecond 
photoionization experiment. The potential energy curves of 
the X, A and E electronic states of Li2 are shown as a func- 
tion of the internuclear distance R as black, green and blue 
solid lines respectively. The potential curve associated with 
the ground electronic state of LiJ is shown as a red solid line. 
Their associated molecular rotational quantum numbers are 
denoted by Nx, Na, Ne and N+ , while M and A/+ denote 
their projection on the electric field polarization axis. 



or TT. If Aip = the two components are in phase, 
and this corresponds arbitrarily to the storage of the in- 
teger n = 1. If Aip = TT the two components are out of 
phase, and this is interpreted as the storage of the integer 
n = m. 

With two vibrational levels ve, one can store in a sim- 
ilar way the combinations 00, 01, 10 and 11, correspond- 
ing to the integers rt = 0, 1, 2, and 3 where the first binary 
digit is determined by the phase shift between the rota- 
tional components of the first vibrational level while the 
second digit is controlled similarly by the second vibra- 
tional level. This can be generalized, and with Uy vibra- 
tional levels one can store all integers up to n = 2"" — 1. 
Encoding this type of classical information in a quantum 
system is of course a very challenging task, and efficient 
schemes are necessary to reveal this information in the 
last experimental step. We propose here to read this in- 
formation with the help of a picosecond ionizing pulse, 
as shown in Figure ^ 

This last laser pulse, the probe pulse, has a bandwidth 
Ao; which encompasses the two rotational components 
associated with the various vibrational levels without 
mixing the different vibrational states during the ioniza- 



tion process, i.e. 

2B,ot{2NA + l) < hAuj < hu^ib , (1) 

where Brot and Wvib are the molecular rotational constant 
and the vibrational frequency in the i?-state. Section Hvl 
will show that the encoded number can be revealed by 
an interference effect modifying the time-resolved photo- 
electron spectrum P{£) which exhibits a series of peaks 
characteristic of the recorded integer. 

III. THEORETICAL APPROACH 

A. The molecular basis set 

We follow the dynamics of the Li2 molecule by prop- 
agating in time its associated electro-nuclear wave func- 
tion \E'(re , i?, i) decomposed in two parts 

*(ri,^,t) = *B(ri,^,t) + «'+(ri,ii,t) (2) 

corresponding to the i5-electronic state of Li2 and to 
(Li J -I- e~) respectively. The coordinates of all electrons 
are denoted by the vector rj, , and the vector R = (R, R ) 
represents the internuclear coordinate. We now sepa- 
rate the global electronic coordinate t|. of all electrons 
into the coordinate rj. of the core electrons and the co- 
ordinate f of the ionized electron. We then write the 
electro-nuclear wave functions '5 e and as the follow- 
ing Born-Oppenheimer expressions 

^E{re,R,t) = yjEiR,t)^E{re\R) (3a) 
*+(ri,il,i) = yj+{r,R,t)<^>+{r,\R) (3b) 

where $_E(?e |^) and $+(rc \R) denote the electronic wave 
functions associated with the i?-state of Li2 and with the 
ground electronic state of Li^ respectively. 

The i?-electronic state (^S+ symmetry) is now con- 
sidered as a 3scr Rydberg state, and the electronic wave 
function <i>£;(re|i?) is expressed in the molecular frame 
(Hund's case (b) representation) as 

'i>E{r,\R) = Mr,ra\R)Yooir) (4) 

On the other hand, the electronic part of ip+i'^', R,t) is 
expressed in the laboratory frame (Hund's case (d) rep- 
resentation) as 

^+ {r,R,t)^JdeY, t) M^,r\R) y,,„ (f ) (5) 

£.771 

where (j}£{£,r\R) is the electronic continuum wave func- 
tion at energy e and {£, m) denote the electron angular 
momentum and its projection in the laboratory frame. 
Note that even though the one-photon transition consid- 
ered here from the i^-electronic state will result in the 
ejection of a p-electron with i = 1, we keep the double 
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sum over £ and m in the following for the sake of gener- 
ality. Various similar approaches have already been used 
for the calculation of femtosecond time-resolved photo- 
electron angular distributions T7|. 

In the present study the photoionization step is per- 
formed by a picosecond laser pulse of linear polarization. 
The following selection rule therefore applies for the pro- 
jections M and Af+ of the molecular rotational angular 
momenta of L\2 and LiJ on the polarization axis 



M = K'V 



(6) 



For any given initial projection Af , an unambiguous re- 
lation therefore relates M"*" and m. Hence the quantum 
number Af + is replaced in the following by M — m. 

In our time-dependent approach a propagation is per- 
formed independently for each initial value of Af, and the 
angular parts of the nuclear wave packets ipEiRjt) and 
ip'l^{e, R,t) are thus expressed as the following expan- 



(7a) 



Ne 



N+ 



in terms of the normalized Wigner rotation matrices 



V^j j^{R) which verify 



where the nuclear Hamiltonian TLfj is given by 



2n 



52 NsiNE + l) 



i?2 



Ve{R) (12) 



We assume that the i5-state coherent superposition of 
rovibrational levels has been created with a femtosecond 
laser pulse which has been sent thought a phase mask 
such that the phases fvE,NE associated with each level 
can be controlled. The binary information is then stored 
in the phase differences Aipy^^ = ^ve-Na+i ~ Vve,Na-i 
between the two rotational components Ne = Na ± 1 of 
each vibrational level ve- In addition we have also as- 
sumed that an amplitude mask has been used to com- 
pensate the Franck-Condon factors which normally gov- 
ern the A X transition. The initial populations of the 
different levels in Eq. (|9a|l have therefore been taken as 
independent of the vibrational quantum number ve- 

The molecule is then submitted to an ionizing laser 
pulse associated with the linearly polarized electric field 



E{t) = Eq f{t)cos{ujt)e 



(13) 



where Eq and lo denote the electric field amplitude and 
the angular frequency of the radiation, e is the unit po- 
larization vector. The pulse envelope f{t) is defined by 
the Gaussian-like expression 



dR V^'f, {R ) Pf/, ^f,, [R ) = 5nn' 5mm- Saa' (8) 
where 6 stands for the Kronecker delta symbol [T^ . 



B. The laser- molecule interaction 

The initial electro-nuclear wave packet (O is prepared 
at time t = from a "source level" {vx , Nx) as a coherent 
superposition of several rovibrational levels (ve,Ne) in 
the two-step process depicted in Figure ^ The initial 
components of the expansions iQ are thus given by 

<-(i?,0) = cneM J2 e'^""'"" XvE,NEiR) (9a) 

Ve 

</_:!(£, i?,0) = (9b) 
where 



CNe,M — V^Ne + 1 {2Na + 1) V^Nx + 1 



Na 1 Nx 
M -A/ 



Ne 1 
A/ 



Na 



Na 1 Nx 


Ne 1 Na 




(10) 



The rovibrational eigenstates Xve.Ne{R) of total en- 
ergy E{ve,Ne) are solution of the time-independent 
Schrodinger equation 



fit) 



(14) 



where r is the Full Width at Half Maximum (FWHM), 
and 2t the total pulse duration. 

Introducing the expansions ©-© and in the time- 
dependent Schrodinger equation, and projecting onto the 
electronic and rotational basis functions yields, in the 
dipole approximation, the following set of coupled differ- 
ential equations for the nuclear wave packets V'm^ [Ri t) 
and ip'^/_:^{e,R,t) 



E{t) ^Y. ^n"eM ^M-L (15a) 

N+,i,m 

E{t)y,Y.^N"EM' <^ (15b) 

Ne 



where the nuclear Hamiltonian is given by 

^2 N+[N++l) 



a. 

2iJ, 



dR^ 



i?2 



+ V+{R) (16) 



The matrix elements A4 



N+,. 
Ne.M 



which couple the nuclear 



'^Ne Xve-Ne{R) — E{ve,Ne) Xve.Ne{R) (11) wave packets evolving on the electronic potential curves 
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Ve{R) and V+{R) read 



M 



Ne,M 



N 



Ne 1 N 
M -M 



N+ £ N 
M — m m — Af 



{Ne\N\N+£} (17) 



where the total angular momentum (ion + electron) 



7V= TV 



(18) 



has been introduced. 

The source term {Ne\N\N^£} is evaluated in the 
molecular frame. The electron wave function originally 
described in the laboratory frame in Eq. © is expressed 
in the molecular frame (Hund's case (b) representation) 
using the frame transformation technique implemented 
by Ugo Fano in his pioneering work on Ha pj|. Follow- 
ing these lines, we obtain 

1 

{Ne\N\N+£} = {N+£\NA) e"^^ 



A=0 



X di{e,R) {NeI\NA) (19) 



where A is the projection of the total angular momentum 
N on the molecular axis. 

From the definition given in Eq. IjlSj) . the following 
relation holds 



A = A+ + A 



(20) 



where A+ and A are the projections of the ion and elec- 
tron angular momenta N and £ on the molecular axis. 
In the present case, the ionic core presents a sym- 
metry (A+ — 0) and hence 



A = A 



(21) 



In Eq. H19|l . the phases (tt/ia) with A = or 1 represent 
the phase shifts of the a and tt electron continuum wave 
functions relative to the regular radial Coulomb function. 
The short range quantum defects /ia with A = or 1 are 
indeed associated with the S and 11 p-Rydberg series of 
Li2 [13. 

The matrix element {N^£\NA) of the unitary frame 
transformation and the Honl-London rotational factor 
{Ne1\NA) of Eq. lO are defined by 



{N'£'\NA) = (-1) 



TV' 



' N N'' 



-A A 



(22) 



In the sum of Eq. ((T^ d^{e,R) denotes the energy and 
i?-dependent ionization dipole moment from the £^-state 



4(£,i?) 



£ 1 
-A A 



di{e,R) 



(23) 



where 

di{s,R) = J </>!(£, r|i?)$;(r,|i?) 

X r (j)E{r,fc\R) drdfc (24) 

The present study is limited to a restricted range of 
photoelectron energies (e < 200 cm~^) and of internu- 
clear distances (in the vicinity of the _B-state equilibrium 
distance Re)- This justifies the Condon approximation 
di{e,R) ~ est used hereafter. The quantum defects /is 
and /in are also taken as independent of R, and their nu- 
merical values /is ~ 0.001 and /in — —0.287 have been 
extracted at R = Re from the electronic potential ener- 
gies given in pol |. 



C. The time propagation 

To calculate the ionization of Lia subjected to a pulsed 
laser radiation, we propagate the nuclear wave pack- 
ets 'ip^f{R,t) and 'ilj^^^^{e, R,t) in time during the en- 
tire pulse using the split operator method developed by 
Feit et al HH 



/ 



,N+.e 



\ 



I 



_ ~iHSt/h 



I 



\ 



VM-m 



t+St 



\ 



(25) 



where the total (molecular -I- interaction) Hamiltonian 
Ti = T + V + W{t) is split in three parts corresponding 
to the kinetic (T), potential {V) and interaction {W{t)) 
propagators 

^-iHSt/h _ ^-itst/2h ^-ivst/2h ^~iW(t)St/h 

^^-^vst/2n^-^fst/2n^^^^^3^ (26) 

The kinetic and potential propagations are performed in 
the momentum and coordinate spaces respectively. Fast 
Fourier Transformation (FFT) allows rapid passage back 
and forth from one representation to the other at each 
time step. The propagator associated with the laser in- 
teraction term W{t) is calculated using a simple diago- 
nalization of its associated interaction matrix |22|. 

Since we are only dealing with bound vibrational states 
in this study, typical grids extend from R — 3.5 au to 
R — 14.0 au with 2^ grid points. The potential energy 
curves Ve{R) and V+{R) are taken from ^20]. Numeri- 
cally, the rotating wave approximation (RWA), very ac- 
curate for the present case of low laser intensities and 
vertical resonant transitions, allows a substantial gain 
of computational time. This approximation, which con- 
sists in neglecting the so-called counter-rotating terms, 
results in a simple one-photon dressing (energy transla- 
tion by huj) of the ion electronic potential V-^{R) in our 
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time-dependent approach |23| ■ A time step of the order 
of 5t ~ 4 fs ^ T is then sufficient for convergence. 



D. The photoelectron spectra 

The analysis of the photoelectron angular distribu- 
tions is made by projecting the wave packet -0+ de- 
fined in Eq. jSJ at the end of the pulse {t — 2t) on 
the energy-normalized solutions of the field-free ionized 
molecular states l^^. We need to define the set of out- 
going plane waves elastically scattered in the direction 
k= {k,k) for a prescribed asymptotic electron kinetic 
energy e = h^k"^ /2m. These are represented by the usual 
expansions on angular momentum states |24| 



thus giving 



\e.k) 



Y;^{k)Me,r\R)Yernir) (27) 



where £,e denotes the Coulomb phase shift of each partial 
wave. We therefore evaluate, for a given value of M, the 
angular distribution of the ejected photoelectron at some 
fixed energy e as 



PM{e,k)= dRdk {e,k\ilj+{f,R,2T))f 



(28) 



This multiple integral can be written in a more conve- 
nient form with the help of the expansions Q and ((7b)). 
After the integration over the electronic coordinate 
f= {r,r) and over the angle R = {9r, 0^) we get 



PM(£,fc) =^ ^ ^z(^'~^)e^«-«^')r;„(fc)r,„.(fc) 



X jdR (e,i?,2T) Vf/ji(e,i?,2r) (29) 

One can notice here the appearance of an incoherent sum 
over the quantum numbers and m, while the sum 
over the electron angular momentum € is coherent and 
gives rise to various {I x €') cross-terms. The total pho- 
toelectron angular distribution at some fixed energy e is 
then obtained by averaging over the initial distribution 
of M 



P(e,fc) (x^PM(e,fc] 



(30) 



M 



In the present study, a p-electron is ejected and the 
value of the electron angular momentum is therefore fixed 
io t =1. In this case, Eqs. (^5)) and simphfy to 



p{e,k) 



M N+ rn 



v. 



Af+,1 
M—m\ 



(e,i?,2T) dR (31) 



The total photoelectron spectrum is then obtained by a 
summation over the ejection angle k = {Ok, 4>k) 



(32) 



P(e) = / dkP{e,k) 



M N+ rn 



dR (33) 



Numerically, the continuous variable e is discretized in 
150 energy values, with 10 ^ e < 190 cm^^. The proba- 
bility that the electron exits in the k direction with an 
energy e is calculated from Equation (|31|l on a grid of 
325 points in fc = (Okjipk)- 

Even if Eq. (|33(l only reveals a series of incoherent 
sums, one should not forget that each exit channel N'^ 
may be reached from different initial levels Ne- As a con- 
sequence, an initial quantum superposition of rotational 
levels can induce an interference effect in the photoelec- 
tron spectra which arises from the phases of the different 
components in this initial wave packet. 



IV. RESULTS AND DISCUSSION 

A. Analysis of the photoelectron spectra and of 
the interference effect 

Typical photoelectron spectra are shown Figure|21when 
a single initial rovibrational level {ve,Ne) is populated. 
The upper and lower graphs of Figure |21 correspond to 
the same initial vibrational state ve ^ 0, but to dif- 
ferent initial rotational excitations: Ne = 1 in the up- 
per panels (a) and (b), A^^; = 3 in the lower panels (c) 
and (d). It is assumed here that these levels are pre- 
pared using the two-step process depicted Figure ^ with 
{vx = 0,Nx ^ 1) and {va ^0,Na = 2). The incoher- 
ent sum over M in Eqs. H31|l and (|33|) therefore extends 
from M = — ltoM = l only. The photoelectron spec- 
tra shown in the left panels (a) and (c) have been cal- 
culated using the accurate short range quantum defects 
/is = 0.001 and iijj = —0.287 of Li2 while the right pan- 
els (b) and (d) correspond to an hypothetically isotropic 
Li2 molecule with jis = fJ-n- 

The solid blue lines in Figure [3 represent the pho- 
toelectron spectra obtained with the laser frequency 
u! = 14184 cm~^ (wavelength 705 nm) and with the pulse 
duration r = 15 ps. One can notice that for the initial 
condition Ne = 1 (panel (a)) two peaks are seen in the 
photoelectron signal. Their energies correspond to the 
ion exit channels = with iV+ = 1 or 3. 

The V-^. — ve vibrational selection rule is observed in 
this case because of the pronounced Rydberg character 
of the i?-electronic state at short internuclear distances. 
The two potential curves Ve{R) and V+{R) are indeed 
parallel for i? ^ 7 au (see Figure . Around i? ~ 9 au, 
an avoided crossing gives rise to the double-well struc- 
ture of the i?-state, deeply modifying the Rydberg na- 
ture of this excited state. However, this avoided crossing 
has a negligible influence on the lowest vibrational levels 
{VE ^ 4). 
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FIG. 2: (Color online) Photoelectron spectra calculated us- 
ing a single initial state {ve,Ne) as a function of energy e. 
Upper panels (a) and (b): the initial state is ve ~ 0, Ne = 1. 
Lower panels (c) and (d): the initial state is ve = 0, A''^; = 3. 
The photoelectron spectra calculated using the quantum de- 
fects of Li2 are shown in the left panels (a) and (c), while the 
spectra in the right panels (b) and (d) correspond to /ie = ^J.^ 
(see text for details). The solid blue lines correspond to the 
pulse duration (FWHM) r = 15 ps while r = 2.5 ps for the 
red dashed lines. The laser wavelength is 705 nm. The photo- 
electron energies expected from a simple energy conservation 
rule e = E{ve, Ne) + Huj — E{v+, N^) with the usual selec- 
tion rule v+ = Ve are shown as thin vertical dotted blue lines 
for various ion rotational quantum numbers A''^ . 



The branching ratio between the iV+ = 1 and ~ 3 
exit channels clearly favors iV+ = Ne — 1. This result 
is not unexpected since in the case of an isotropic sym- 
metry (corresponding here to /is = /in) the other ioniza- 
tion channel {N~^ = 3) has no intensity, as one can see 
in the panel (b) of Figure |2| The molecular anisotropy, 
expressed by the phase difference (/is — /in)7r — 0.37r, al- 
lows for the exchange of angular momentum between the 
ionized electron and the nuclear rotation. This effect, 
which can be seen as a mutual electron-nuclei entangle- 
ment, explains the appearance of the additional branch 
N~^ = 3 in panel (a). 

For the initial level Ne — 3 (lower part of Figure [21, 
this effect is also seen with the appearance of the two 
satellite peaks N~^ = Ne ± 2 around e = 45 cm"-'^ and 
around e = 60 cm"-'^. These two peaks are not seen in 
panel (d) when /^s = /^n- Note also from panels (a) 
and (c) that the central peak, which corresponds to 
an ionization without exchange of angular momentum 
(iV"*" = Ne), is located at the same photoelectron energy 
e — 55 cm~^ for both initial rotational levels Ne = 1 
and 3. This happens because the rotational constants 
of the i?-state of Li2 and of the ground electronic state 
of LiJ are almost identical. 

The dashed red lines shown Figure |21 finally represent 
the photospectra calculated with the same parameters 



except for a much shorter pulse duration r = 2.5 ps. 
In this case the peaks associated with different values 
of the ion rotational quantum number TV"*" overlap due 
to the large spectral bandwidth of the pulse. A single 
broad photoelectron peak is, therefore, obtained around 
e = 55 cm^^ whatever the initial rotational level Ne- 

In the isotropic case shown on the right hand side of 
Figure |21 the photoelectron peaks calculated with this 
shorter pulse duration have a symmetric shape since a 
single ionization channel (A''"'" = Ne) is observed. On the 
other hand, a slightly asymmetric shape is obtained with 
the real Li2 molecule (see the red dashed line in panel 
(c) for instance), due to the asymmetric distribution of 
the two satellite peaks = Ne ± 2 on both sides of the 
central peak corresponding to N~^ = Ne ■ 

In the experiment, it is expected that two Ne rota- 
tional levels can be populated in a coherent distribu- 
tion whose initial phase difference is controlled using the 
phase mask of a pulse shaper. This type of experiment 
has been implemented recently for higher rotational lev- 
els in the Group of Stephen Leone 13] for instance. The 
photoelectron spectra calculated using in-phase and out- 
of-phase initial distributions of Ne = 1 and 3 are shown 
in FigureOlas red solid and blue dashed lines for the pulse 
duration r = 2.5 ps. One can recognize in these two spec- 
tra the slightly asymmetric distributions discussed pre- 
viously. 



(a) T = 2.5 ps 



— A<p = 

- - A<p = It 

A<p = it;tij, = fijj 



IS 



(b)T= 15 ps 




e (cm ) 



FIG. 3: (Color online) Photoelectron spectra calculated us- 
ing a single initial vibrational state {ve ~ 0) with a coherent 
superposition of two rotational levels (A''^; — 1 and Ne ~ 3) 
as a function of energy e. The pulse duration is r = 2.5 ps 
in the main graph (a) and t = 15 ps in the small upper-right 
inset (b). The laser wavelength is 705 nm. The red solid line 
and blue dashed line correspond to an in phase {Aip = 0) and 
an out of phase {A(p = n) coherent superposition respectively. 
The isotropic case /ie ~ /in is represented by the green dotted 
line for Aip — and the green solid circles for Aip = n. The 
labels TV"*" ^ Ne in the inset assign the various energy peaks 
with respect to the initial (Ne) and final (A'^"'') rotational 
levels. See text for details. 

A clear and significant interference effect is also seen 
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here. The total ionization probability is amplified by a 
factor of 2.05 when comparing the in-phase and out-of- 
phase initial conditions. 

Because of the orthogonality of the Wigner rotation 
matrices jy{R) the photoelectron probability can be 
written as an incoherent sum over the rotational quantum 
number 7V+, as one can see in Eq. The interference 

effect seen here is therefore not due to an interference 
between the different exit channels labeled by iV+ , but is 
due to an interference taking place in the same ionization 
channel between the two pathways coming from the two 
possible initial levels Ne- 

Let us consider for instance the ionization channel as- 
sociated with = 3. This final state can be reached 
through the main branch coming from Ne = 3 or through 
the small satellite arising from Ne = 1. One could think 
at first sight that the cross term associated with this in- 
terference mechanism is probably negligible considering 
the small branching ratio (~ 0.1) between the N^ = Ne 
and N^ — Ne ± 2 pathways (see Figure However, 
even in this apparently unfavorable case, it can be easily 
estimated that this cross term can induce a significant 
interference effect which could in principle modify the 
ionization probability by a factor of three when compar- 
ing destructive and constructive interferences. 

In our case all exit channels are not subjected to this 
interference. For example N'^ — 5 can be reached from 
Ne = 3 only. In addition the photoelectron signal is a 
complex incoherent average over the different possible 
values of M and m (see Eq The two different path- 

ways Ne N+ and N'^ N+ are also affected by the 

different values of the coupling matrix elements -^^^'^™ 
given in Eq. (|17|l . The contrast separating a constructive 
from a destructive interference is therefore not maximum, 
and after this complex averaging Figure |31 shows that the 
interference effect changes the ionization probability of 
Li2 by a factor of about two. 

Our analysis is also confirmed by the photoelectron 
spectrum obtained with /is = /in (green dotted line and 
green solid circles in Figure OJ. Whatever the phase shift 
A(p between the Ne = 1 and Ne = 3 initial rotational 
levels the same photoelectron spectrum is obtained, and 
no interference effect is seen. This is due to the disap- 
pearance of the Ne N^ — Ne ± 2 pathways, as dis- 
cussed previously. 

The interference effect can also be suppressed by in- 
creasing the pulse duration r. In this last case the photo- 
electron peaks associated with different rotational quan- 
tum number N~^ do not overlap anymore, and two identi- 
cal N^ originating from different Ne never appear at the 
same energy (see the small inset in Figure|31for instance). 
As can be seen in Eq. (|33|) . the wave packets associated 
with different photoelectron energies never interfere and 
as a consequence the interference effect seen with shorter 
pulses is not observed with long pulses. We have indeed 
verified that the spectrum shown in the small inset of 
Figure|ni(T — 15 ps) does not depend on Aip. 

We will now show that the interference effect shown 
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FIG. 4: (Color online) Photoelectron spectra calculated using 
two initial vibrational states {ve ~ and 1) with a coherent 
superposition of two rotational levels (A''^; — 1 and Ne ~ 3) 
as a function of energy e. The pulse duration is r = 2.5 ps 
and the laser wavelength is 699.8 nm. The dashed blue 
lines show the photoelectron spectra P{n, e) calculated when 
the integer values n = (panel a), 1 (panel b), 2 (panel c) 
and 3 (panel d) are stored in the phase difference between 
the Ne = Na ± 1 rotational components of the initial wave 
packet. The solid red lines in the small insets show the signal 
difference S{n, e) = P{n, e) — P(0, e) for n = ... 3. 



in Figure 13 with relatively short pulses can be used to 
reveal efficiently a more complex phase information ini- 
tially stored in the rotational wave packet of Li2 in a way 
similar to the implementation performed in the group of 
Philip Bucksbaum with Rydberg wave packets [lOl iH^ . 



B. Applications to quantum information 

In Li2 , we use the quantum superposition of two ro- 
tational levels in different vibrational states to store a 
binary information. The case of two vibrational levels 
is shown in Figure^ The combinations 00, 01, 10 and 
11, associated with the decimal values n — 0, 1, 2 and 
3, can then be stored. The measurement of this infor- 
mation is performed by photoionization, thanks to the 
anharmonicity of the electronic potential curves of the 
molecule. 

Since the vibrational frequency of the i?-state and of 
the ion are almost identical, in an harmonic approxima- 
tion and with = ve the energy of the ionized electron 



e = E{vE,NE) + nuj- E{v+,N+) 

harm 



(34a) 



est + Srot [NEiNE + l) - N+{N+ + 1)] (34b) 
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does not depend on the value of the initial vibrational 
level. However, the anharmonicity of Li2 is large enough 
(~ 20 — 50 cm^^) to exceed by far the rotational spacing 
(Slot — 0.5 cm^^). Different initial vibrational levels can 
therefore be clearly distinguished even with the relatively 
short pulses (1 ps ^ r ^ 5 ps) which allow for the obser- 
vation of the interference effect discussed previously. 

The dashed blue lines in Figure ^ represent the pho- 
toelectron spectra calculated with the phase differences 
(0,0), (0,7r), (7r,0) and (7r,7r) in the vibrational states 
{ve — ^,ve ^ 0) between the rotational levels Ne = 1 
and Ne — 3. These phase differences correspond to the 
integers n = 0, 1, 2 and 3. We denote by P{n,e) the 
associated photoelectron spectra. The interference effect 
(amplification of the ionization probability by a factor of 
two) is clearly seen for both vibrational levels. 

If we define the signal difference S{n, e) by the follow- 
ing relation 



Sin,e) = P{n,e) - P{0,e) 



(35) 



a direct visualization of the binary representation of n 
is obtained in the graphs showing S{n,e) as a function 
of the energy e. These signal differences are represented 
as red solid lines in the four small insets of Figure 01 
This photoelectron signal difference S{n, e) character- 
izes the cross terms associated with the — Ne and 
7V+ = N'^ ± 2 paths. It can also be seen as a direct mea- 
sure of the electron-nuclei entanglement induced by the 
molecular anisotropy. 

In order to check the scalability of this method, we 
have calculated the signal difference S{n,e) for n = 0, I, 
10 and 31 when the binary information is stored in the 
first five vibrational levels of the i^-state. These results 
are shown as red solid lines in the four panels of Figure 
One can notice the effectiveness of the method for such 
a small number of vibrational levels. 

A few tendencies are worth being noted in this figure. 
Since the total (electron+ion) energy is fixed, the lowest 
energy peaks can be assigned to the highest ion vibra- 
tional levels w+. The fact that the energy separation 
between two peaks increases with is simply explained 
by the increasing anharmonicity of the potential curve 
with larger internuclear distances. 

Finally, one can also notice that the energy peak on 
the left is not as intense as the other ones. This peak 
corresponds to an ion left in t;_|_ = 4 and therefore to the 
molecule being initially in ve = 4. This initial level is al- 
ready located at an energy close to the avoided crossing of 
the i?-state (see Figurenj, and the v+ = ve selection rule 
is no longer fully verified. This initial level is effectively 
also ionized in — 1 = 3, yielding the ejection of 

an electron carrying much more energy. This effect ex- 
plains the smallest probability observed in Figure for 
the Ve = 4: ^ v-^- — A pathway. 

Higher vibrational levels also show this tendency of 
being ionized into vibrational levels with < ve- This 
type of behavior could be seen as a limitation of the pro- 
posed mechanism for storing and reading information in 
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FIG. 5: (Color online) Photoelectron spectra calculated us- 
ing five initial vibrational states {ve = ... 4) with a coherent 
superposition of two rotational levels {Ne ~ 1 and A*'_b — 3) 
as a function of energy e. The pulse duration is r = 2.5 ps 
and the laser wavelength is 699.6 nm. The solid red lines 
show the signal difference S{n, e) — P{n, e) — P{0, e) when 
the integer values n — (panel a) , 1 (panel b), 10 (panel 
c) and 31 (panel d) are stored in the phase difference be- 
tween the Ne = Na ± 1 rotational components of the initial 
wave packet. The step function shown as a green solid line 
in panel (d) represents the vibrational Franck-Condon factors 
^(e) = {v+\ve) associated with the different photoelectron 
peaks. The ratio S{n, e)I ,^(e) is shown in panel (d) as a blue 
dashed-hne. 



and from the rotational degree of freedom of a diatomic 
molecule. However, this effect can easily be corrected 
for by renormalizing the photoelectron signal with the 
vibrational Franck-Condon factors 



xUR) XvAR) dR 



(36) 



between the eigenfunctions associated with the ve and 
vj^ = Ve vibrational quantum numbers. The step-like 
function ^(e) is shown as a green solid line in panel (d) of 
Figure[Sl The decreasing probability seen in the left peak 
associated with = = 4 is compensated by calculat- 
ing the simple ratio S(n, e)/ ^{e). This ratio is shown in 
the same graph as a blue dashed line, and all construc- 
tive interference peaks now reach the same height. Since 
the vibrational wave functions XvJyR) XvJ\R) can be 
easily calculated from the potential curves given in |20j |. 
this correcting procedure is quite straitforward. 
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V. CONCLUSION 

We have proposed here a theoretical model for the 
study of the short-pulse photoionization of Li2 , follow- 
ing a scheme relatively close to recent experimental im- 
plementations. This model takes into account the molec- 
ular rotational degree of freedom and allows for the cal- 
culation of the photoelectron spectrum resolved in time, 
energy and in angle. 

We have used this time-dependent model to predict a 
new and efficient mechanism for measuring binary clas- 
sical information initially stored in the phases of a rota- 
tional wave packet. A picosecond laser excitation induces 
an interference between these different rotational compo- 
nents. This typically molecular interference effect mea- 
sures the electron-nuclei entanglement which takes place 
thanks to the anisotropy of the diatomic molecule, and 
transfers the rotational phase information in the photo- 



electron spectrum. 

These results indicate that a high degree of control can 
be achieved in this type of molecular systems using sim- 
ple pulse shaping techniques. A systematic exploration 
of the control achievable in the photoelectron angular 
distributions will be presented in another paper. 
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